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' We investigate how additive wealc noise (correlated as well as uncorrelated) modifies the param- 

, eters of the Gray-Scott (GS) reaction diffusion system by performing numerical simulations and 

(N . applying a Renormalization Group (RG) analysis in the neighborhood of the spatial scale where 

biochemical reactions take place. One can obtain the same sequence of spatial-temporal patterns 
by means of two equivalent routes: (i) by increasing only the noise intensity and keeping all other 
model parameters fixed, or (ii) keeping the noise fixed, and adjusting certain model parameters to 
' their running scale-dependent values as predicted by the RG. This explicit demonstration validates 

the dynamic RG transformation for finite scales in a two-dimensional stochastic model and provides 
' < I ' . further physical insight into the coarse-graining analysis proposed by this scheme. Through several 

D ' study cases we explore the role of noise and its temporal correlation in self-organization and propose 

I a way to drive the system into a new desired state in a controlled way. 
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INTRODUCTION. 



O , The dynamics of certain biological systems frequently follows some self-organization process where the development 
of new, complex structures takes place primarily in and through the system itself. This self-organization is normally 
' triggered by internal variation processes, which are usually called "fluctuations" or "noise" , that have a positive 
^ : influence on the system. For instance, recent theoretical studies and experiments with cultured glial cells and the 
0^ Belousov-Zhabotinsky reaction have shown that noise may play a constructive role on the dynamical behavior of 
spatially extended systems H], 1^, Therefore, one can not ignore the role of noise in chemical and biological 
self-organization and its relationship with the environmental selection of emergent patterns Q . 

In these systems there is a strong interplay between the equations of the microscopic dynamics and the environmental 
fluctuations acting on the smallest scales. In order to get a qualitative idea of the system behavior one can average out 
the influence of noise on the shortest scales by increasing, infinitesimally, the scale of observation. This can be done 
■ applying dynamic Renormalization Group (RG) techniques 0, . One obtains then the so-called RG flow equations 
' which are generally used to study the RG flow of the system parameters and to find the fixed points in the long 
^ wavelength limit (at the fixed points the system is invariant under the scaling transformation). But our interest in 
I this paper is in the finite-size renormalization (at finite wavelengths or finite scales), i.e. on the infiuence of noise at 
' ^ , the scale where the chemical or biological reactions take place. 

Our final aim is to validate the RG procedure as a general tool which can be used at the scales where these 
biochemical reactions take place and predict the transitions induced by noise with different correlation properties. 
As an example we will investigate the influence of noise on a reaction diffusion system in the vicinity of a transition 
.5^ ' point. Performing numerical simulations and applying the RG analysis at flnite scale ranges we can get a better 
understanding of how weak noise modifies the parameters of this dynamical system and the patterns it converges to. 
^ We consider the Gray-Scott (GS) reaction diffusion model 0, which is one of the simplest models of biochemical 

5^ , relevance leading to spatial and temporal patterns when diffusion is included. Numerical simulations of this system 
in the deterministic case have revealed a surprisingly large set of complex and irregular patterns as a function of the 
system parameters 8|. In addition we want to include the influence of the fiuctuating environment. Fluctuations 
of internal (thermal) origin represent microscopic degrees of freedom, and since they evolve in spatial and temporal 
scales much shorter than those of the gross variables of the system they are assumed to be additive and uncorrelated 
in space and time. When the origin of the noise is external there may exist a coupling between the system and the 
fluctuations. In this case there is no difference between the time and length scales of the noise and the field, and this 
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noise could be multiplicative and have some structure in space or time. We study the following stochastic system 
with additive noise: 

dV 

— = XUV^ - fiV + D.V^V + Tj4x,y,t) (1) 
dU 

— - ua-XUV^-i^U + D^V^U + r,^ix,y,t) 
at 

where = + -g^. We will consider Gaussian white zero- mean noises (uncorrelated) 

< Vu,v{x,y,t)r]u,v{x' ,y' ,t') >— 2Au^v5{x — x')5{'y — y')5{t — t') and Ornstein-Uhlenbeck (OU) correlated noise 

< ilu,vix,y,t)r]u^vix' ,y' ,t') >= Au^v exp ^ S{x — x')S{y — y') (the Gaussian white noise case limit 2Au^v5{t — t') 
is obtained when the correlation time r ^ while keeping = A„^„t constant). In the absence of noise, Eq. ^ 
defines the Gray-Scott model, which is a variant of the autocatalytic Selkov model of glycolisis, corresponding to the 
following chemical reactions: 

U + 2V ^ iV (2) 
V ^ P 
U ^ Q 



V{x,y,t) and U{x,y,t) represent the concentrations of the chemical species U and V, and are functions of the two- 
dimensional space and time variables. A is the constant rate of the reactions, P and Q are inert products, ^ is the 
decay rate of V and ;/ the decay rate of U. The equilibrium concentration of U is uq/v, where is the feed rate 
constant. The chemical species U and V can diffuse with independent diffusion constants and D^. All the model 
parameters, including the noise amplitudes A^^vi are positive. 

In the case of Gaussian white noise the stochastic system described in Eq. ^has been numerically integrated for a 
particular set of parameters Q . The pattern to which the system converged changed drastically with small changes 
in the noise intensity, suggesting that this system can suffer noise induced transitions when subjected to random 
fluctuations [l0|. 

Additionally, one can apply the RG transformation to the Gray-Scott model subject to additive random fluctuations 
The dynamic RG procedure consists of two steps: an elimination or thinning out of the fast or short wave-length 
modes (A < < A where s = e' > 1 is the scale in wavenumber space above which the fluctuations have been 
integrated over, A = ^ and L the minimal length where the equations under consideration are valid) followed by 

a rescaling of the remaining modes (fc — > sk) 0. Applying this infinitesimal transformation one gets the effective 
(i.e. renormalizcd) parameters of the system when going through an infinitesimal scale change. This analysis leads 
one to consider a set of coupled differential equations governing the RG fiow in parameter space when this system is 
influenced by noise. At one-loop order in the noise amplitude, only two out of the total of eight model parameters 
run with scale, namely, the decay rate v, of the U fleld and the nonlinear coupling A, and the corrections are driven 
by the w-noise (not the u-noise). In particular for the case of white Gaussian noise, in a two-dimensional space, the 
pertinent equations read 

dv \AyK2K'^ 
dX / AXA„K,,A^ 



and for the case of OU noise with long-range correlations read: 
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zv + CSC 



(1 + 20„)^ 



(Ai + i?.A2)i+2e. 



(4) 



{2x + z)X 



-4 CSC 



A 



with Ay ~ AyT and 9^ — 1 Here dl is an infinitesimal change in the scale s of the system (fi-om s = e' to 

s = s + ds with ds = sdl). The scaling parameters z and x ^-re the so called dynamic exponent (for time) and 
roughness exponent (for the fields U and V) and are not needed for the present study. 

These equations ()3I4|I predict how the effective parameters of the system run with scale and noise amplitude and 
have been derived and analyzed in detail for the long- wavelength limit (i.e. in the large scale limit s oo) in [lll |. 
Our aim in this paper is to validate these equations in the intermediate scale range, far from the fixed points, where 
the chemical processes take place. To carry this out, we will use H3I4|I to calculate the effective parameters of the 
GS system for a finite and noise of significant intensity. We then solve ^ using these effective parameters (and 
negligible noise), and compare the resultant patterns based on these running parameters with those that are obtained 
numerically when all the model parameters are held fixed and only the noise intensity is increased. 



II. DYNAMIC RG THEORY WITH FINITE SCALING AND UNCORRELATED NOISE. 

The RG consists of thinning the degrees of freedom followed by a re-scaling of length and time 5^ . It is a general 
calculational scheme for treating problems where fluctuations at many length scales are important. The final aim 
of the RG is to describe how the dynamics of a system subject to random fluctuations evolves as one changes the 
scale of observation. In practice, the RG procedure permits one to establish rigorous correspondences between sets 
of parameters defining physically different states observed at different scales. Indeed, (|3I4|I express the fact that v 
and A change with the scale I, as functions of the noise and the other model parameters, as a consequence of the 
coarse-graining. In particular, in the case under consideration, this correspondence is due to the presence of noise and 
a finite change in the spatial scale. 

Following Eq. ^ we observe, that for a scale variation of Al the predicted running or modification of the i/ and 

A parameters induced by noise alone, should be Aiy = +Aljj^^^^ and AA = —Al A (^ (^^.j a'2)\'^+1+d' a'^+d a^) ) 
respectively. For the two-dimensional case, K2 — l/(27r). In our numerical simulation the smallest scale of observation 
is that of the cell with Ax x Ay and Aa; = Ay. This determines the minimal wavelength and maximal momentum 
A = 1^ that can be considered in the problem. Let us take the cell as our initial scale, our scaling factor is therefore 
s = 1. The diffusion term allows us to average data over neighboring cells, this in turn defines a greater scale of 
approximately 3Aa; x 3Ay which leads us to s = s -I- As with As — 2 — sAl and A^ — 2. 

We will explore the patterns that are formed within a certain range of the parameter space. Similar studies have 
been performed in the deterministic Gray-Scott model, keeping A = 1 constant, varying /i between 0.06 and 0.14 and 
v between 0.01 and 0.07 (in this region there is a transition from two stable steady states to one trivial state and a 
great variety of patterns) . 

In Table we show, for a given noise intensity Ay , the predicted effective parameters A -I- AA and i' + Av based 
on the one-loop RG calculations in for \ = \^Dy — 0.5, = I, uq = 0.05 = v and ^ = 0.1155. The labels in 
column 1 and 5 will be used later to name the study cases. In Table UTI we show the same calculations, in a different 
region of parameter space, when the unperturbed values are v = uq = 0.03 and fi = 0.086. In both cases we have set 

= Uq so that the trivial state solution is always U = I. 

The numerical simulations of system evolution have been performed using forward Euler integration of the finite- 
difference equations following discretization of space and time in the stochastic partial differential equations. The 
spatial mesh consists of a lattice of 256 x 256 cells of size Ax = Ay = 2.2, with periodic boundary conditions. Noise 
has been discretized as well. The initial conditions consisted of one localized square pulse with {U — 0.5, V — 0.25) 
perturbing the trivial steady state (U — 1,V = 0) plus random Gaussian noise. The perturbing pulse measured 
22 x 22, just wide enough to allow the autocatalytic reaction to be locally self-sustaining. The system has been 
numerically integrated for up to 30000 time steps (step size At = 1). In the figures, only the concentration of the 
substrate U is shown. When displayed in color, the blue represents a concentration between 0.2 and 0.4, where the 
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TABLE I: Effective or renormalized parameters (sixth and seventii columns) induced by a wfiite noise of intensity wfien tfie 
unperturbed or bare values are fi — 0.1i55, v = uo = 0.05 and A = f . 



case 




A 


V 


case 


A + AA 




A 


4.8 X lO"*" 


f 


0.05 


a 


1 (-fo-'') 


0.05 (+3 X fO"") 


B 


f.6f X 10-" 


f 


0.05 


b 


0.99997 


0.050i 


C 


4.36 X fO"'' 


f 


0.05 


c 


0.9999f 


0.05027 


D 


4.84 X fO"'' 


f 


0.05 


d 


0.999903 


0.0503 



TABLE II: Effective or renormalized parameters (sixth and seventh columns) induced by a white noise of intensity when 
the unperturbed or bare values are = 0.086, f = uo = 0.03 and A = f . 



case 




A 


V 


case 


A + AA 


V + I\v 


I 


4.8 X 10"" 


1 


0.03 


i 


1 (-10-'^) 


0.03 (+3 X 10"") 


II 


1.60 X lO"'' 


1 


0.03 


ii 


0.99997 


0.0301 


III 


4.33 X 10~* 


1 


0.03 


iii 


0.999914 


0.03026 



substrate is being depleted by the autocatalytic production of V ^ yellow represents an intermediate concentration of 
roughly 0.8 and red represents the trivial steady state (JJ — 1). 

A. Noise Induced Transition 

We will now solve numerically the time evolution of this system Q holding all parameters fixed and study the 
resulting patterns as a function of increasing noise intensity A^. First, for A = 1,151, = 0.5, = \, — 0.05, /i = 
0.1155, V = 0.05 and — 0, we integrate this system with very weak noise A^ = 4.8 x 10^^ and obtain the pattern 
shown in Fig. ^A. Then we increase only the noise intensity with the strength given in the second column of Tabled 
(A and v are held fixed to the values given in the next two columns) and obtain the patterns shown in Figs. to 
D. The system evolves first forming long stripes. As we vary only the noise, some of the stripes split into spots and 
finally spots are the only stable structures. Noise induces a transition between different patterns or states. 

Next we repeat this study at a different location in parameter space: uq = v ^ 0.03, fi — 0.086 keeping A = 1, Dy — 
0.5, Du = 1 and setting the noise intensity A^ to the values given in the second column of TablellTI (A and i/ are held 
fixed to the values given in the next two columns). We obtain now the patterns shown in Figs, to III. For weak 
noise, the system fills in the space with a pattern of short tubular structures and many spheres in clusters. Then, as 
we increase the noise intensity, the spheres clusters tend to disappear and the tubular structures become longer. 

We have chosen two particular examples where weak noise induces an observable transition between two patterns (in 
some other cases weak noise or the equivalent small parameter variation does not modify significatively the resulting 
pattern with respect to the reference case). 

B. Transitions induced by parameter variation. 

Next we apply a fixed weak noise of intensity Ay — 4.8 x 10^^ and numerically solve the system ^ with the 
renormalized parameters, A = A + AA and i> = v + Av, that correspond to the coarse-grain scale Al — 2 and the noise 
intensities studied above which are given in the sixth and seventh column of Tables HI and ITU In Fig. [21-(a) to (d), we 
show the results for uo = 0.05 and ^ = 0.1155 and in Fig. 01-(i) to (ui) for uq = 0.03 and fi — 0.086. The results are 
equivalent to the patterns obtained before when only the noise was applied. Within the limits of perturbation theory, 
one can therefore obtain the same sequence of patterns and morphologies either by increasing the noise intensity alone 
or by holding the noise fixed and, applying a coarse-graining analysis, adjusting the A and parameters to their 
renormalized values as predicted by the RG. This is the main result of this paper. 

Note that there is some discrepancy between the study cases (d)-(D) and (iii)-(III): in case (d), where only the 
deterministic parameters are changed, no self-replicating structure survives, whereas in (D), with noise, the pattern 
persists and in case (iii) the spheres have almost disappeared whereas in (III) the structures are erased. This difference 
may indicate that we are reaching the threshold where the higher order corrections (beyond one-loop) of the RG can 




FIG. 2: Results obtained with renormalized parameters A + AA and v + Ai^ as given in Table |T] and holding the weak noise 
fixed = 4.8 X lO"*^. Figures (a)-(d) are shown at t = 30000. 



no longer be ignored, or the point where the GS system is not renormalizable (new effective reactions are generated 
which can not be described by a simple renormalization of the original parameters as was explained in 11]). We have 
observed that generally, for « 4.5 x lO"** and beyond, the patterns can become too noisy and/ or non equivalent 
to the patterns based on the renormalized parameter calculation. The dimensionless perturbation parameter of the 
RG analysis performed in is 5 = [i^+d^k-^) which when evaluated at this limit gives g « 4.5 x 10^''. 

Although the equivalence given by the RG equations breaks down at a certain noise intensity, the change that a 
pattern will suffer when applying weak noise (for instance, from long stripes into spots or from many spheres into 
long tubes) is easily predicted. 



III. DYNAMIC RG THEORY WITH FINITE SCALING AND CORRELATED NOISE. 

Next, we include a study case of this system with temporally correlated noise. Following Eq. |0J) we observe, that for 
a scale variation of A/ the predicted modification of the v and A parameters induced by noise alone, should be Ai/ = 
+^U,,il±2l^_^0^ and AA = -A^csc \A^K,K\u + D^K^) (;;^^;;r;:|::g:^:;:^'"" A, re- 

spectively (where 0^ = 1). 

Our starting point now is equivalent to the previous study case (c), A — 0.99991 and v — 0.05027, which (with 
the weakest white Gaussian noise A^ — 4.3 x 10^"*) produces spots, see Fig|Sl-(a). Next we add some OU noise with 
Av = AyT = 2.61 X 10~^ and r = 2000, see FigEK/?), under this condition the system produces stripes and spots. 
Increasing the OU noise intensity further to A^ — A^t — 7.67 x 10^"^ and r — 5876, see Fig|31-(7y), the pattern shows 
only long stripes. Adding OU noise of increasing strength to a system operated with the parameters of (c) we have 
reproduced the sequence (c) to (a). 

We can explain this behaviour as follows: this correlated noise has induced a variation in A and v of opposite sign 
to the case of white Gaussian noise leaving the system with new effective parameters A = 1.0001 and = 0.05 for the 
study case (77). If we solve the system with these renormalized parameters we recover the same situation as in the (a) 
study case which produces only stripes. For the intermediate case (/3) the OU noise has induced a variation to the 
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FIG. 4: Results obtained with renormalized parameters A + AA and v + as given in Table UTI and with a fixed weak noise 
= 4.8 X 10"''. Figures (i)-(iii) are shown at t = 20000. 

new effective parameters A = 0.99991 and v = 0.05027. We have included for completeness a simulation with these 
values (and the weakest white Gaussian noise = 4.3 x 10~^), see Fig|Sl-(7), the resulting pattern with stripes and 
spots is equivalent to the (/3) case as expected. 

Next we choose as starting point the study case (iii), v — 0.03026 and A = 0.999914, see FigEl-(l), which produces 
long tubes and few spheres. We add some OU noise with A„ = A^t = 5.22 x 10""^ and t — 2000, see Fig|SK2), 
under this condition the system produces more spheres and shorter tubes. This correlated noise induces a variation 
in the effective system parameters: A = 1.00001 and D — 0.030072. If we now solve the system with the renormalized 
parameters we obtain again a pattern with more spheres and shorter tubes as in the study case (2), see Fig|Sl-(3). 

Applying correlated noise one can also find noise induced transitions. The RG prediction for the parameter variation 
is also valid for the case of correlated noise. Note that noise does not always drive the system in the same pattern- 
sequence direction since the sign of parameter variation/renormalization depends on the type of noise. 

IV. CONCLUSION 

We have studied the dynamics of the GS system with uncorrelated and correlated noise. The environmental 
fluctuations, the spatial diffusion and the non-linearities of the system create cross-coupling terms and the new 
effective system parameters depend on the noise amplitude and correlation time. We have shown that for weak 
noises the lowest order one-loop RG analysis can be applied at intermediate scale ranges, and be used to estimate the 
effective parameters of the system. In other words, by combining analytic and numerical work, we have established an 
equivalence between a sequence of patterns generated by varying the noise amplitude but keeping all other parameters 
fixed and a companion sequence generated by keeping the noise fixed and varying (i.e., renormalizing) instead some 
of the model parameters according to the RG flow equations. 

For suitable noise intensities the perturbed values give rise to new patterns. Noise does not always drive the system 
in the same pattern-sequence direction since the sign of parameter variation depends on the type of noise. Based on 
this knowledge and using weak noise, one can lead the system into a new state in a controlled way. Notice also that 
adding noise to a system which is expected to produce structures of a certain spatial correlation length (such as for 
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FIG. 5; Results for correlated OU noise, (a) is the starting point with fi = 0.1155, uo = 0.05, A = 0.99991 and u = 0.05027 (/?) 
is obtained adding intermediate OU noise, (77) is obtained with stronger OU noise. Adding OU noise is equivalent to operating 
the system with new effective parameters: for the intermediate study case {(3) this is shown in (7) which has been solved for 
A = 0.99996 and = 0.05018, for case [i]) the induced parameters are A — 1.0001 and (> — 0.05 as in the previously shown 
study case (a). The patterns are given at t = 30000. 

(1) (2) (3) 




FIG. 6: Resuhs for correlated OU noise. (1) is the reference case for fi = 0.086, uo = 0.03, A = 0.999914 and u = 0.03026 
(2) is obtained adding OU noise and (3) changing the system parameters to the new effective parameters A = 1.00001 and 
v = 0.030072. The patterns are shown at t = 20000. 



instance the spots in the study case (a)) does not necessary drive the system into the production of patterns with 
shorter spatial correlations. On the contrary, it can drive the dynamics of the system into a new situation which 
produces patterns with longer spatial correlations (such as the stripes in the study case (77) with OU noise). 

This raises important questions on the role of noise in chemical and biological self-organization and the environ- 
mental selection of emergent patterns. When describing a biological system not only should we estimate the relevant 
parameters of the system, but also the magnitude of the stochastic influences, the spatial scales involved and the 
correlation function. The interdependence of most biological systems with the environment makes this result specially 
relevant. In the case of non-linear systems and in the vicinity of transition points, the role of noise can be non-trivial, 
allowing the system to explore new, well organized situations which might be favored by the external natural selection 
processes. Applying the dynamic RG theory at finite scales one can have an intuitive idea of how the system will be 
driven by the environmental fluctuations. 
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